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Abstract 

In the framework of the LSDA+[7 method electronic structure and magnetic properties of the in- 
termetallic compound Gd2Fei7 for both rhombohedral and hexagonal phases have been calculated. 
On top of that, ab initio exchange interaction parameters within the Fe sublattice for all present 
nearest and next nearest Fe ions have been obtained. It was found that for the first coordination 
sphere exchange interaction is ferromagnetic. For the second coordination sphere exchange inter- 
action is observed to be weaker and of antiferromagnetic type. Employing the theoretical values 
of exchange parameters Curie temperatures Tc for the first as well as for the first and second 
coordination spheres of both hexagonal and rhombohedral phases of Gd2Fei7 within Weiss mean- 
field theory were estimated. Obtained values of Tc and its increase going from the hexagonal to 
rhombohedral crystal structure of Gd2Fei7 agree well with experiment. Also for both structures 
LSDA+f/ computed values of total magnetic moment agree well with experimental ones. 

PACS numbers: 75.50.Ww, 78.20.-e, 71.20.-b 
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I. INTRODUCTION 



Intermetallic compounds i? 2 Fei 7 with high concentration of iron (R is a rare-earth ion) 
are of interest for both experimental and theoretical investigations.- These compounds are 
distinct for a relatively large magnetic moment - mainly caused by the magnetic moments at 
the Fe ions - and comparatively low Curie temperature Tc- The compounds are crystallized 
in two related crystal structures - rhombohedral T/i 2 Zn 17 -type and hexagonal T/i 2 iVz 17 -type. 
In each of these structures, the Fe ions can be located in four inequivalent crystallographic 
positions, and for this reason various experiments reveal four values of the local magnetic 
moments at different positions of the Fe ions. Also, most of the intermetallic compounds 
-R 2 Fei 7 have collinear magnetic structure except for Ce 2 Fei 7 , Tm 2 Fei 7 , and Lu 2 Fei 7 ,- where 
the magnetic moments of the Fe ions in all crystallographic positions are collinear to each 
other, but frequently antiparallel to the magnetic moments of the rare-earth ions R. 

Initially, the moderate Curie temperatures Tc in the compounds i? 2 Fei 7 was explained 
on the base of assumption 3 of a critical distance r c « 2.5 A between two ions in a Fe-Fe pair, 
so that for the distances r < r c the exchange interaction is antiferromagnetic, whereas for 
r > r c - ferromagnetic. From this assumption it followed that a low value of Tq is caused 
by a competition of ferromagnetic and antiferromagnetic exchange interactions between 
neighboring pairs of the Fe ions in various crystallographic positions. It also meant that an 
expansion of the lattice should result in an increased number of ferromagnetic interactions 
and corresponding growth of Tc- Such an expansion of the lattice can be reached either 
inserting nonmagnetic interstitial atoms,-^ or substituting the Fe ions by nonmagnetic atoms 
with the larger ionic radii.- 1 ^ 

However, in a number of experiments, for example, for substitutional solid solutions 
-R 2 Fei 7 _a;Si x , it was found that the crystal lattice contracts, and Tc grows .-^^ These results 
question the adequacy of the simple model of exchange interactions of neighboring exchange- 
bounded ions that change the sign at the critical distance r c between the ions. 

To determine a relation between Tc and exchange interactions in the _R 2 Fei 7 series, the 
compound Gd 2 Fei 7 is of especial interest. First, among the series Gd 2 Fei 7 has the highest 
Curie temperature Tc- Secondly, Gd itself is located at the middle of the rare-earth se- 
ries, and Gd 2 Fei 7 separates the compounds with light rare-earth elements crystallizing in 
the rhombohedral T/i 2 Zn 17 -type structure, from the compounds with heavy rare-earth ele- 
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ments crystallizing in the hexagonal Th2Nin-type structure. Thus Gd2Fei7 is an intriguing 
example among the i?2Fei7 series which can be found in both crystal structures. 

Apparently, in most experiments samples are mixture of two structures. At the same 
time, there are experimental works for either pure rhombohedral structure^ 11 or pure 
hexagonal structure^ of Gd2Fei7. For Gd2Fei7 in the rhombohedral structure, the val- 
ues of critical temperature and saturation magnetization were found to be T^ h =475 K,— 
Mg h =21.5 /xe/f.u.,— respectively. For the hexagonal structure, these values are T^ ex =466 K 
and Mg ex =21.2 /ig/f.u.— Although these two sets of experimental data differ slightly and 
further experimental study is needed, one can notice that Gd2Fei7 in the rhombohedral 
phase has the higher values of magnetic characteristics, than in the hexagonal one. 

Motivated by this fact, in the paper we consider how the transfer from the rhombohedral 
to hexagonal structures in Gd2Fei7 influences the values of the local magnetic moments and 
exchange interaction parameters. The later values are obtained using the results of elec- 
tronic structure calculations within LSDA+[7 method. Further values of Curie temperature 
were estimated in the framework of the Weiss mean-field theory. Moreover based on these 
calculations we propose explanation of rather low Tc values and relatively small saturation 
magnetization of i^Fe^ series with respect to elemental bcc Fe. 

II. CRYSTAL STRUCTURES AND LOCAL MAGNETIC MOMENTS 

The intermetallic compound Gd2Fei7 can crystallize both in the rhombohedral ThiZnyi- 
type structure (space group R3m - no. 166 in International Tables for Crystallography) 
and the hexagonal Th%N z^-type crystal structure (space group P63/HIHIC - no. 194). Main 
structural blocks are iron hexagons building layers. Most of these hexagons have empty 
centers, whereas others contain the rare-earth Gd ions or mediate the Fe ions in dumbbell 
positions. 

The rhombohedral structure of the Th 2 Zni 7 -tjpe is shown in Fig. [TJ In this type of 
structure the Fe(18/) hexagons (let us denote Fe(18/) as Fe3 for the rhombohedral structure) 
contain the Gd ions in positions 6c either at the center of each hexagon, as it is presented in 
the upper and lower layers in Fig. (TJ or alternate with empty hexagons, as it is shown in the 
inner two layers. Interlayer Fe(6c)=Fel ions at the dumbbell positions lower the symmetry 
of some iron ions in the nearest intermediate layers from 18h to 9d (above and below the 
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dumbbell), and hence, these intermediate layers without Gd contain ions Fe(9cf)=Fe2, as 
well as Fe(18/i)=Fe4 in the hexagons. Unit cell of the rhombohedral structure consists of 
one formula unit of Gd2Fei7. 

The hexagonal structure of the T/i 2 ^i7-type is shown in Fig. [2J In contrast with the 
rhombohedral structure, there are no layers where each hexagon contains the Gd ion. In 
this structure the Gd ions are located only at the centers of the hexagons (let us denote 
Fe(12j)=Fe3 for the hexagonal structure) alternating with the hexagons with the dumbbells 
of Fe(4/) (Fe(4/)=Fel), and the layers with the hexagons containing Gd at each center are 
absent. Moreover, here one can find the alternating intermediate layers Fe(6g)=Fe2 and 
Fe(12fc)=Fe4 without the Gd ions. Unit cell of this structure consists of two formula units 
of Gd2Fei7 due to the different symmetry of the rare-earth ions coordination environment 
in the aft-plane translations. 

For both structures in the layers with alternating Gd-containing and empty hexagons Fe3, 
these hexagons are deformed - the hexagons with Gd are slightly expanded, the hexagons 
without Gd - contracted. Furthermore, because for the hexagons Fe2 and Fe4 the interme- 
diate layers the Gd ion and the dumbbell Fel ions are alternating from one layer to another, 
so the ions Fe2 and Fe4 in these layers are slightly distorted upward and downward with 
respect to the plane, i.e. these layers are slightly corrugated. 

The electronic structure of the rhombohedral structure of Gd2Fei7 was calculated in the 
LSDA+t/ method (Ref. Il3l ) in the framework of the band calculation package TB-LMTO- 
ASA (tight-binding, linear muffin-tin orbitals, atomic sphere approximation). 14 However let 
us mention here that only Gd 4f-shell was treated within the LSDA+[7 method since cor- 
relation effects are much stronger there than for the Fe 3d-shell (Ucd ~7 eV (see below), 
Up e ~2 eV— ^). We believe that greater part (not all) of electronic interaction effects is 
taken into consideration for Fe 3d-shell in the frame of LSDA approximation. It is well 
enough justified since responsible Stoner parameter for Fe 3d-shell Jp e ~1 eV is close to 
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Up e . For the hexagonal structure, the results of such a calculation were reported in Ref. 
below we include some details of that calculation to provide the full picture for compar- 
ison of the rhombohedral and hexagonal structures of Gd2Fei7. Lattice parameters for 
the rhombohedral structure are a=8. 538 A and c=12.431 A (Ref. Q), for the hexagonal 
structure - a=8.496 A and c=8.341 A (Ref. Il2l ). Atomic spheres radii were chosen as 
i?(Gd)=3.72 a.u. and _R(Fe)=2.66 a.u. for the hexagonal structure, and i?(Gd)=2.86 a.u.— 



4 



and i?(Fe)=2.62 a.u. - for the rhombohedral structure.— Orbital basis contained 6s, 6p, 5d, 
and 4/ muffin-tin orbitals for Gd and 4s, 4p, and 3d for Fe sites. Integration over the first 
Brillouin zone was performed using 32 irreducible k-points (6x6x6=216 is a total number 
of k-points). Empty atomic spheres without nuclear charge were inserted in the case of the 
hexagonal structure to fill the interstitial regions. 

In the LSDA+£7 method Coulomb interaction for Gd 4f-orbitals was taken into account 
via parameters of direct U and exchange J Coulomb interactions of the 4 /-electrons in Gd. 
For Gd2Fei7 in both structures the calculations of these parameters in the constrained LDA 
method^ resulted in £7^=6.7 eV and Jcd=0.7 eV that is in agreement with the values for 
Gd metal.^ 

Local magnetic moments for various Fe and Gd sites obtained in the calculations are 
listed in Table U for both phases under consideration. Note that in qualitative discussions 
of magnetic properties of i?2Fei7 especial attention is frequently paid to the iron ions in the 
dumbbell positions (Fel in our notations). From Table [J one can see that magnetic mo- 
ments of the Fel ions for the hexagonal structure are larger than for the rhombohedral one. 
However, total magnetic moment M(SFe) value of the iron subsystem in the rhombohe- 
dral phase is larger (M r / l (EFe) =38.04 /ig/f.u. vs. M/j ea ,(SFe)=36.54 fi B /f.u.. Experimental 
value of total Fe ions magnetic moment for the mixture of two phases is 36.9 /xs/f.u.).— 
Magnetic moments of the Gd ions in both calculations are antiparallel to the iron magnetic 
moments, so that the total magnetic moments are M^EFe + Gd)=23.70 /ig/f.u. and 
Mfe ex (SFe + Gd)=22.21 /ig/f.u., respectively. These values agree well with experimental 
ones Mg h —21.5 /i B /f.u. n and Mg ex =21.2 /i B /f.u., 12 correspondingly. 

III. EXCHANGE INTERACTIONS 

In the framework of spin-fluctuation theory of magnetism it was founcP^ 1 ^ 1 ^ that the 
magnetic behavior of transition metals can be described by the Heisenberg model with long- 
range exchange interactions between classical spin vectors, in general, their magnitudes are 
non-integer. For an intermetallic magnet with several types of magnetic moments a spin 
Hamiltonian of such a model can be written as 
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TABLE I: Calculated values of local magnetic moments for both rhombohedral and hexagonal 



phases of Gd2Fei7. 



rhombohedral structure 


hexagonal structure 


one 


A/T ( a 
M [flB) 


oirc 




Gd(6c) 


-7.17 


Gd(26) 


-7.13 






Gd(2d) 


-7.20 


Fe(6c)=Fel 


2.19 


Fe(4/)=Fel 


2.31 


Fe(9d)=Fe2 


2.26 


Fe(6#)=Fe2 


2.10 


Fe(18/)=Fe3 


2.17 


Fe(12j)=Fe3 


2.40 


Fe(18/i)=Fe4 


2.31 


Fe(12A;)=Fe4 


1.87 



H= -\H E ^(la-n 6 )S a (l a )S 6 (n 6 ), (1) 

a,b l a ,n b ,\ a ^n b 

where a and b are types of magnetic ions (e.g. a=Fel, Fe2, Fe3, and Fe4), l a - radius 
vector of a magnetic ion of the a type, S a (l a ) - classical spin vector with the magnitude S a , 
corresponding to magnetic moment fi a = g^sSa, and I a b(\a — n £>) is an exchange parameter 
between ions of the types a and b with the distance l a — between the ions. 

Later an ab-initio method to calculate the exchange interaction parameters between dif- 
ferent lattice sites for the classical Heisenberg model at T=0 was proposed.— In that scheme 
the exchange interaction parameters are determined calculating the second-order derivative 
of the total energy at T=0 with respect to small deviations of magnetic moments of the 
corresponding lattice sites from collinear magnetic configuration. This method was used in 
the present work to calculate the exchange interaction parameters I a b (a, 6=Fel, Fe2, Fe3, 
Fe4) in the iron sublattice of both structures of Gd2Fei7. 



First note that in Ref. 
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the energy of exchange interaction between different 
sites of the magnetic lattice is presented in the form E ex (l — n) = — 2J^(l a — 
n b )e a (\ a )e b (n b ) as an exchange between two classical unit spin vectors e a (l a ) = 
(siri(j) a (l a )sinQ a (l a ), cos(p a (l a )sinQ a (l a ), cosQ a (l a )) and e b (n b ) regardless of the magnitude 
of the magnetic moments fi a and fi b at the sites l a and rij, under consideration. Since the 
classical spin vector S a (l a ) is equal to S a (l a ) = S a e a (l a ), hence the relation between the ex- 
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change parameters J^ b calculated using the method of Ref. [24j and the exchange parameters 
I a b of the Hamiltonian (Tj[|) will be given as 2 J^ b = I a bS a Sb. Moreover, some authors^ calcu- 
lating exchange parameters for insulators prefer to write the energy of exchange interaction 
as E ex {\ — n) = — J^f(l — nt)S a (l a )S5(n fe ) via the exchange parameters between classical 
spin vectors with the magnitude S a , Sb=l/2. It is easy to demonstrate that if one uses the 



method of Ref. [25j to calculate J^f , then our exchange parameters I a b are related with 
via the equality I ab = J$ /(4S a S b ). 

Generally, in the intermetallic compounds of the R,2Fen type there are many types of 
exchange interaction: R-R, R-Fe, and Fe-Fe. The weakest interaction is an indirect inter- 
action of the R-R type between the rare-earth ions, because even for the nearest rare-earth 
ions their wave functions practically do not overlap. 

Exchange interaction .R-Fe is also indirect since direct overlap of 4f orbitals of the R ions 
and 3d orbitals of the neighboring iron ion is absent, and only small overlap of a polarized 
5d orbital of the rare-earth ion and 3d orbital of the iron ions is possible. Theoretical 
estimation of this interaction is still difficult, nevertheless, from some indirect experiments 
this interaction can be estimated as 10-30 meV.— The sign of the exchange interaction R-Fe 
defines the orientation of the magnetic moment on the rare-earth ion with respect to the 
magnetic moments of the iron sublattice that seriously affects the saturation moment. At 
the same time, the dominating contribution to free energy of the magnetic system is due to 
the interaction of the Fe-Fe type which defines in fact the value of the critical temperature 
Tc- Because the number of iron ions is 8.5 times greater than the number of R ions, and in 
the Fe subsystem - there is also direct exchange between the neighboring Fe ions and the 
overlapping 3d wave functions. 

In such complicated crystall structures as Th^Znxi and ThzNin the nearest neighbors in 
different crystallographic directions are arranged at different distances. However, analyzing 
the nearest neighbors of magnetic ions of each type Fel, Fe2, Fe3, and Fe4 in both structures, 
the following regularities can be found. If one considers a spherical layer with the radii 
from r[ h =2.385 A < r < ^=2.740 A around a particular ion Fea (a=l, 2, 3, 4) in the 
rhombohedral structure, then all nearest magnetic neighbors will be within this coordination 
layer. Next is an empty spherical layer without ions ^=2.740 A < r < r^ h =4.044 A. The 
second spherical coordination layer will be at 4.044 A < r r 3 h < 4.3 A. 27 The only exception 
from this regularity is the Fe3 class with the nearest neighbor Fe3-type ion at r[^=3.563 A 
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between the first and second coordination layers. It is typical that the atomic sphere of 
the chosen ("central") ion Fea (a—1, 2, 3, 4) and atomic spheres of nearest Fea (a=l, 2, 3, 
4) ions overlap. Therefore the direct exchange gives a dominant contribution in exchange 
parameters of the nearest neighbor Fe ions. 

Similarly, in the hexagonal structure of Gd2Fei7 each ion Fea (a=l, 2, 3, 4) can be 
surrounded by the first spherical coordination layer with the radii rj te:r =2.399 A < r < 
rg =2.735 A, then an empty spherical layer without ions with ^=2.735 A < r < 
r hex = 4 iQg anc [ second spherical coordination layer with 4.108 A j r\ ex \ 4.3 A^I. 
In exactly the same way the Fe3 ions violate this regulation, since these ions have a neigh- 
bor at r^ x =3.605 A between the first and second coordination layers. 

The results for the exchange interaction parameters I a bi)-a — ^b) between different pairs 
of iron ions (in units of K) are presented in Tables [TTHVl In Table [III the parameters of 
exchange interaction (in the descending order of absolute values) between the central ion 
Fea (a=l, 2, 3, 4) and ions of different types from the first coordination layer are presented. 
Here 134 (1) stands for the value of the exchange parameter between ions Fe3 and Fe4 at 
the nearest distance between these ions (index 1 in the brackets), then ^34(2) denotes the 
value of the exchange parameter between ions Fe3 and Fe4 at the second-order distance 
between these ions (index 2). Moreover, in Table HT1 the number of neighbors Z^ir) for a 
given distance are presented. This value gives the number of neighbors for the central ion 
Fea at the distance r. One should note that in our crystal structures Z a b(r) 7^ Zf, a (r) (e.g. 
Z i3 (l)=6 and Z 3 i(r)=2, i.e. the Fel ion has 6 neighbors of Fe3 at ri 3 (l)=2.735 A, whereas 
the Fe3 ion has only 2 Fel ions at the same distance r3i(l)=2.735 A). 

It is observed that for the hexagonal Gd2Fei7 all exchange parameters for each Fea ion 
with other ions Fe from the first spherical layer are positive i.e. ferromagnetic. Surprisingly, 
the exchange interaction parameter between the Fe4 ions within the corrugated hexagon 
plane J44(l)=218.5 K is just a little less than the exchange parameter Jn(l)=238.8 K be- 
tween the Fel ions at the dumbbell positions. Also, one should note that there is no simple 
dependence between the distance of the interacting Fe ions and the value of exchange inter- 
action. For example, in the pair Fe3-Fe3 at the second-order distance r 33 (2)=2.512 A the 
exchange interaction parameter is equal to /33(2) = 123.4 K, whereas in the same plane of 
hexagons Fe3 in the pair Fe3-Fe3 at the distance r 33 (2)=2.403 A the exchange interaction 
parameter is smaller, namely /33(1)=80.4 K. 
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TABLE II: Parameters of exchange in the hexagonal structure of Ga^Fe^ for the ions of the first 
coordination sphere. 



N Exchange (K) Distance (A) Number of neighbors 



1 




=238.8 


m(l) 


1=2.400 


*n(l)=l 




2 




=218.5 




1=2.477 


Z44(l)=2 




3 




=136.1 




1=2.511 


2 3 4(1)=^4 3 (1)= 


=2 


4 




=123.4 


r 3 3(2) 


1=2.512 


«33(2)=1 




5 


/24(1) = 


=101.9 


r24(l) 


1=2.459 


«24(1)=4, 242(1) 


i=2 


6 


/ 23 (i) 


=90.7 


r 23 (i; 


1=2.464 


z 23 (i)=4, z 32 (i; 


i=2 


7 


^13(1) 


=90.5 




1=2.735 


z 13 (i)=6, z 3 i(i; 


i=2 


8 


/ 33 (1) 


=80.4 


r 33 (i; 


1=2.403 


^ 33 (i)=i 




9 


W) 


=72.0 




i=2.653 


Zl 4 (l)=3, Z4l(l) 


i=l 


10 


A 2 (l) 


=69.3 




1=2.616 


«12(1)=3, «21 (1) 


i=2 


11 


Im (2) 


=61.3 




1=2.672 


z 34 (2)=z 43 (2)= 


--2 



Type 

Fel (dumbbell) - Fel (dumbbell) 
Fe4 (corrugated plane) - Fe4 (corrugated plane) 
Fe3 (dumbbell) - Fe4 (dumbbell) 
Fe3 - Fe3 
Fe2 - Fe4 (corrugated layer) 
Fe2 (corrugated layer) - Fe3 (upper and lower layer) 
Fel (dumbbell) or Fe3 - Fel (dumbbell above and below) 

Fe3 - Fe3 

Fel (dumbbell) - Fe4 (nearest corrugated layer) 
Fel (dumbbell) - Fe2 (nearest corrugated layer) 
Fe3 - Fe4 (corrugated layer above and below) 



The method of exchange interaction parameters calculation employed here works in such 
a way that does not distinguish different mechanisms. There are several exchange interac- 
tion scenario super exchange, double exchange, direct exchange and RKKY one. First two 
require nonmagnetic interstitial ions in the system, that is not the case here. Thus only 
last two mechanisms are applicable. Direct exchange mostly come from orbitals overlap. 
One can suppose that absence of straightforward correlations between exchange parameter 
values and corresponding interatomic distances for the first coordination sphere for different 
crystallographic directions is connected with the way how Fe-3d orbitals overlap in these 
directions. However classification of the orbitals for these quite complicated crystal struc- 
tures is highly non-trivial task which we postpone for future investigations. In addition 
some anisotropy in the first coordination sphere of exchange parameters might come from 
RKKY interaction since the Fermi surface of these compounds is rather complicated. For 
the second coordination sphere RKKY interaction is dominant but still might be highly 
anisotropic because of direction dependence of kp. 

In Table ILTT1 the parameters of exchange interaction of the ion Fea (a=l, 2, 3, 4) with 
some ions from the second-order coordination layer with r > r3 ea: =4.108 A, as well as 



9 



TABLE III: Parameters of indirect exchange in the hexagonal structure of Gd2Fei7 for ions of the 
second coordination sphere. 



N 


Exchang 


;e (K) Distance (A) Number of neighb 


ors 


Type 


1 


^34(3)= 


-14. 


^34(3) 


/I 1 or 

=4.135 


2: 34 (3)=2, Z43(3)= 


--2 


Fe3 - Fe4 


2 


^34(4)= 


=5.b 


^34 (4 j 


A 1 

=4.172 


/ A \ O f A\ 

z 34 (4)=2, 243(4)= 


--2 


Fe3 - Fe4 


3 


I M (2)= 


=4.2 


ri4(2J 


=4.201 


zi 4 (2)=3, 241 (2)= 


-1 


Fel (dumbbell) - Fe4 (corrugated layer) 


4 


I23 (2)= 


=-7.9 


^3(2; 


A AOO 

=4.038 


2 2 3(2)=4, 2 32 (2)= 


--2 


Fe3 - Fe2 (corrugated layer) 


5 


43(3)= 


-10.7 




=3.605 






Fe3 - Fe3 


6 


/ 33 (5) = 


-13.2 


r 33(3) 


=4.257 


2 33 (5)=2 




Fe3 - Fe3 


7 


/ 12 (2)= 


-14.6 


n 2 (2) 


=4.108 


Z12(2)=3, 221 (2) = 


--2 


Fel (dumbbell) - Fe2 (corrugated layer) 


8 


/ 34 (5)= 


-16.3 


^34 (5) 


=4.220 


2 34 (5)=2, 243(5) = 


--2 


Fe3 - Fe4 (corrugated layer) 


9 


/ 22 (2)= 


-18.1 


r 22 (2) 


=4.260 


2 22 (2)=4 




Fe2 (corrugated layer) - Fe2 (corrugated layer) 


10 


/ 33 (4)= 


-29.0 


r 3 s(4) 


=4.239 


233 (4) =2 




Fe3 - Fe3 (layer above and below) 


11 / 22 (1)=- 


-111.1 


r 22 (!) 


=4.182 


z 22 (l)=2 




Fe2 - Fe2 (layer above and below) 



the exchange parameters 133(3) for ^3(3) = r^ x =3.605 A in a Fe3-Fe3 pair at the third- 
order distance, in the hexagonal structure of Gd 2 Fei 7 are presented. At these distances the 
exchange parameters are mostly negative, antiferromagnetic, and are usually an order of 
magnitude less than the exchange parameters with the ions of the first-order coordination 
layer. At the same time, a sufficiently strong antiferromagnetic exchange i2 2 (l)=~Hl-l K 
in both pairs of the nearest neighbors Fe2-Fe2 from the neighboring corrugated planes is 
present. 

In Table HV] the parameters of exchange interaction of the ion Fea (a=l, 2, 3, 4) with 
the ions from the first coordination layer, as well as the exchange parameters ^33 (2) with 
the Fe3 ions located between the first and second layers, in the rhombohedral structure of 
Gd2Fei7 are shown. Transferring from the hexagonal to the more compact rhombohedral 
structure of Gd2Fei7, the distance in a dumbbell Fel-Fel decreases from rff x (l)=2.400 A to 
r[^(l)=2.385 A. Where from the exchange parameter /n(l) increases from ij\ ea; (l)=238.8 K 
to III (1)=287.5 K. However, the second largest interaction ^44(1) in the corrugated hexagon 
plane and the third largest interaction ^34(1) between the layers of Fe3 and Fe4 are 
smaller in the rhombohedral structure of Gd 2 Fei 7 : J 4 h 4 e3; (l)=218.5 K vs. 7 4 t(l) = 182.2 K. 
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TABLE IV: Parameters of exchange in the rhombohedral structure of Gd 2 Fei7 for the ions of the 
first coordination sphere and indirect exchange 133(2). 



N Exchange (K) Distance (A) Number of neighbors 



1 


hi(l)= 


=287.5 


m(i; 


==2.385 


2 


744(1)= 


=182.2 


r 44 (l) 


=2.490 


3 


7 34 (1)= 


=125.9 


r 3 4 (l' 


=2.551 


4 


7 33 (2)= 


=121.9 


^33 (2; 


=3.563 


5 


7 24 (1)= 


=121.0 


r 2 4 (l^ 


=2.448 


6 


7 34 (2)= 


=105.7 


^34 (2) 


=2.613 


7 


7i 4 (l)= 


=88.8 


n 4 (l) 


=2.639 


8 


723(1)= 


=87.1 


r 23(l^ 


=2.423 


9 




=83.6 


ri2(l) 


=2.602 


10 


ha(l)= 


=74.1 


n 3 (i; 


=2.740 


11 


7 33 (1)= 


=-36.5 


^33(1^ 


=2.466 



Type 

Fel (dumbbell) - Fel (dumbbell) 
Fe4 (corrugated plane) - Fe4 (corrugated plane) 
Fe3 - Fe4 (corrugated plane) 
Fe3 - Fe3 

Fe2 (corrugated layer) - Fe4 (corrugated layer) 
Fe3 - Fe4 (corrugated layer) 
Fel (dumbbell) - Fe4 (corrugated layer) 

Fe3 - Fe2 (corrugated layer) 
Fel (dumbbell) - Fe2 (corrugated layer) 
Fel (dumbbell) - Fe3 
Fe3 - Fe3 



z 44 (l)=2 
Z34(l)=24 3 (l)=2 
«33(2)=1 

z 24 (l)=4, z 42 (l)=2 

^3 4 (2)=^43(2)=2 

zu(l)=3, z 4 i(l)=l 
z 23 (l)=4, z 32 (l)=2 

«12(1)=3, Z2 1 (l)=2 

«i 3 (l)=6, z 31 (l)=2 
2Sa(l)=2 



While rj 4 ex (l)=2.478 A is smaller than r$(l)=2.490 A, and rgf(l)=2.511 A smaller than 
r^(l)=2.549 A 7|f : (l)=136.1 K exceeds 7^(1)=125.9 K. Also, in the rhombohedral struc- 
ture one antiferromagnetic exchange (with negative value) appears since I^ X (1)=80A K 
in the hexagonal structure becomes 7g3(l)=-36.5 K in the rhombohedral one, whereas 
r|f (1)=2.403 A will be enhanced to rg(l)=2.466 A. Similar results are obtained for rhom- 



bohedral Y 2 Fei7 m Ref. 
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Finally, Table |V] contains exchange parameters for the second coordination sphere of 
rhombohedral Gd 2 Fei7. As in the case of hexagonal phase most of them are antiferromag- 
netic and considerably weaker in comparison with the first coordination layer. 

From these data it is clear that a decrease of the distance in a pair of interacting ions in 
the first coordination sphere results in an increase of the ferromagnetic exchange parameter. 
But since the distance in the basic structural elements changes in opposite way from the 
hexagonal to rhombohedral structure (the distance decreases in the dumbbells, and mostly 
increases in the hexagons), then it is reasonable to use the Curie point Tc to estimate 
collective effect of these changes in the distances and corresponding exchange parameters 
transferring from one structure of Gd 2 Fei7 to another. Let us use Weiss mean-field theory 
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TABLE V: Parameters of exchange in the rhombohedral structure of Gc^Fe^ for the ions of the 
second coordination sphere. 



N 


Exchang 


;e (K) Distance (A) Number of neighbors 


Type 


1 


Ji 4 (2)= 


-12.6 


ri4(2j 


=4.206 


zi 4 (2)=3 


Fel (dumbbell) - Fe4 (corrugated layer) 


2 


i23(2)= 


= 1.0 


^23(2) 


=4.063 


z 23 (2)=4,Z32(2)=2 


Fe2 (corrugated layer)- Fe3 


3 


-^34(3) = 


=-1.2 


7-34(3) 


=4.132 


z 34 (3)=z 43 (2)=2 


Fe3 - Fe4 (corrugated layer) 


4 


Ji2(2)= 


=-3.6 




=4.095 


zt)(2)=3 zoi (2)=2 


Fel (dumbbell) - Fe2 (corrugated layer) 


5 


I 34 (5)= 


=-5.0 


ru (5) 


=4.190 


z 34 (5)=z 43 (5)=2 


Fe3 - Fe4 (corrugated layer) 


6 


/ia(2)= 


-15.2 


ri3(2) 


=4.237 


zi 3 (2)=6,z 3 i(2)=2 


Fel (dumbbell) - Fe3 


7 


/ 34 (4)= 


-16.4 


r 3 4(4) 


=4.171 


2 34 (4)=z 43 (4)=2 


Fe3 - Fe4 (corrugated layer) 


8 


/33(3)= 


-20.4 


^33(3) 


=4.192 


2 33 (3)=2 


Fe3 - Fe3 


9 


I 44 (2)= 


-36.7 


r 44 (2) 


=3.841 


z M (2)=l 


Fe4(corrugated layer) - Fe4 (corrugated layer) 


10 


144(3)= 


-57.5 


r 44 (3) 


=4.236 


z 44 (4)=2 


Fe4(corrugated layer) - Fe4 (corrugated layer) 



taking into account only the Fe-Fe exchange interactions. 

Within the Weiss mean-field theory the Hamiltonian (CQ) can be written as 

H MF = ^EE J ^ A ») < S a(M >< *f0. + > " E E (U (2) 

a,b la,Ai, a l a 

where A ab — n& — l a , and molecular field /i a (l a ) for spin S a (l a ) can be presented as 

h a (L) = 7 ^(A a6 ) < Sf (l a + A ab ) >= ^ 7 a6 (| A ab \)z ab (\ A ab \)a b . (3) 

Here cr b =< S b (l a + A ab ) > denotes the thermodynamic average of ^-projection of classical 
spin for the ion Fefe (6=1, 2, 3, 4). z ab (\A ab \) is the number of Fe6 neighbors at the distance 
|A a fc| from the central ion Fea. 

Calculating the average value a a in a field h a , one obtains a set of self-consistent equations 
for a a (a=l, 2, 3, 4): 

n t /haSa \ n r ., /h a S a . kgT 

(L(x) = coth(x) — 1/x - Langevin function). Linearizing the set Eq. (j4]) with respect to 
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small a a , one obtains a set of linear equations for Tc'- 

k B T c a a = yJ2J2 ^(\ A ab\)zab(\ A a6 |K (5) 

b \A ab \ 

Let us solve the set of equations (jSJ) for the rhombohedral and hexagonal structures of 
Gd 2 Fei 7 restricting ourselves by using I ab (\A ab \) and z a b(\A ab \) as the parameters of di- 
rect exchange and number of neighbors in the first coordination layer (Tables HVl and HT1). 
Upon these assumptions the estimations for the Curie point for exchanges only from the 
first coordination layer are T c h =429 K and T^ ex =402 K. These estimations are below the 



experimental values T c h =A75 K (Ref. M) and T£ ex =466 K (Ref. 
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Such a situation is 



typical for estimations of Tq for ferromagnetic metals in the framework of the molecular 
field approximation (c.f. for instance, estimations for Fe, Co and Ni^). At the same time 
these theoretical values grasp a trend of the Curie temperature growth in transferring from 
the r hombo hedral to the hexagonal structures in Gd2Fei7 what takes place in experiments 
Refs. 
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Here we can address also a question how Tc changes in case second coordina- 
tion sphere is taken into account. As shown above in Tables IrlWl next nearest exchanges are 
predominantly AFM-like but are weaker. Thus one should expect a bit of lowering of result- 
ing T c . Indeed, we found T^(lst+2nd)=378 K and T£ e:c (lst+2nd)=353 K. Nevertheless a 
trend of the Curie temperature growth from the rhombohedral to the hexagonal structure 
is still observed. 

But comparing these values one should keep in mind several reasons why the calculated 
results of Tc cannot precisely fit the experiment. 

Firstly, in accord with spin-fluctuation theory of magnetism approximating magnetic ex- 
citations in metals with the classical Heisenberg Hamiltonian, all exchange interactions in 
metals are long-range, relatively slowly decaying, and sign oscillating (similar to Ruderman- 
Kittel interaction).—^ This conclusion also follows from our results for the hexagonal struc- 
ture of Gd 2 Fei 7 , see Table [III] (for the rhombohedral structure see Table fVjl . where the 
exchange interaction of the Fe ions with the neighbors of the second-order coordination 
layer mostly change their sign and become negative, antiferromagnetic. Exact account of 
these oscillating interactions requires either exact functional form of the exchange distance 
dependence for a particular pair of ions (e.g., Fel-Fe3), or calculation of the exchange pa- 
rameters for such a pair on the scale of a few oscillations. At the moment the method of 



Ref. 



24 allows one to calculate the exchange parameters, but it is not clear how to present an 
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explicit functional dependence of the exchange on distance. Nevertheless, in principle in this 
method one can calculate exchange parameters on the scale of a few oscillations, but for a 
metallic magnet with several types of magnetic moments and complicated crystal structure 
it is a cumbersome problem. 

Secondly, also according to spin-fluctuation theory in metals,— the exchange parameters 
of the classical Heisenberg Hamiltonian in the general case are temperature dependent. For 
this reason, the exchange parameters at T=0 and T=Tq differ, whereas in our calculations 
we estimated Tq from the exchange parameters calculated using the T=0 method of Ref. 



24 



Thirdly, an obvious limitation of such calculations is the Weiss mean-field theory, re- 
stricted in its accuracy. 

Before we summarize let us also compare magnetic properties of elemental bcc Fe to ones 
averaged over first coordination sphere of Fe ions of hexagonal Gd2Fei7. First of all value 
of nearest neighbor exchange interaction of bcc Fe (393.2 K, Ref. LL9j; 373.6 K, Ref. 1301) i s 
four times larger than the above mentioned averaged value of Gd2Fei7 (99.6 K, Ref. l3ll ). 
It can be understood from difference of the nearest neighbor distance of bcc Fe (2.482 A) 
and averaged distance of Gd2Fei7 (2.574 A, Ref. 32 ) . Thus hexagonal Gd2Fei7 has more 
loose-packed crystal structure than bcc Fe. This more loose-packed crystal structure also 



leads to a sma 
(1299 G, Ref. 



ler saturation magnetization at T=0 of Fe ion sublattice in hexagonal Gd2Fei7 



34]) in contrast to bcc Fe (1740 G, Ref. [331) . It becomes even more clear since 



values of bcc Fe local magnetic moment (2.2 /zg) is almost equal to one averaged over all 
crystallographic positions for hexagonal Gd2Fei7 (2.15 Similar results can be obtained 
for rhombohedral phase of Gd2Fei7. 



IV. CONCLUSION 



In this work we calculate from first principles electronic structure and values of exchange 
interaction parameters of the Fe ions sublattice for the hexagonal and rhombohedral phases 
of intermetallic compound Gd2Fei7 within first and second coordination spheres. Based 
on the values of exchange interaction paremeters Curie temperateres Tq were also calcu- 
lated. Obtained theoretical Curie temperatures T£ fc =429 K and T^; ex =402 K are found to 
be slightly below the experimental values T^ h =475 K— and T^ ex =466 K— . However the 
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tendency of T c to be higher than Tq is well reproduced. The same is found for the case 
if second coordination sphere is taken into account. Corresponding Curie temperatures are 
T^(lst+2nd)=378 K and T£ ex (lst+2nd)=353 K. 

Beside that from our calculations we observe that for the first coordination sphere ex- 
change between different types of Fe ions is ferromagnetic (with only one exception for 
rhombohedral phase). At the same time exchange for next nearest neighbors is mostly an- 
tiferromagnetic one. Latter one is consistent with RKKY exchange picture (exchange sign 
depends on pair bond length) However such competition only is not enough to 

decrease significantly Curie temperature Tc of i?2Fei7 as proposed if Ref. O. Also it was 
reported before that exchange value between Fe in dumbbell positions is the determinating 
factor for the Tc valued. However in this paper we showed that there are several exchange 
interactions almost of the same strength as dumbbell one, e.g. J44 (see Tab. HT) . 

Let us also roughly analyze values of exchange interactions in terms of crystal struc- 
ture. Collecting data from Tables HIUVl one can observe for both phases analogous tendency. 
Exchanges within any layer are weaker than interlayer ones. Exceptional is only Fe2 ion sur- 
rounded with four Fe4 ions with exchange 124(1). Anisotropy of direct exchange interaction 
parameters in the first coordination sphere comes from Fe-3d orbitals overlap anisotropy. 
On top of that RKKY exchange gives slight modulations in the first coordination sphere 
because of complexity of the Fermi surface of the compounds. For the second coordination 
sphere RKKY is main contribution to exchange interaction but for the reason mentioned 
above is also highly anisotropic. 

After all we propose another one observation to explain relatively low Tc of R2Fei7 series. 
To this end we compare magnetic properties of Fe ion sublattice of hexagonal Gd2Fei7 to 
elemental bcc Fe. There are two ingredients influencing Tc value: local magnetic moment 
values and exchange interaction values. Despite the fact that for both hexagonal Gd2Fei7 and 
bcc Fe local magnetic moments of Fe ion are almost identical, exchange interaction values of 
bcc Fe is four times stronger then for hexagonal Gd2Fei7-~ It comes from more dense packing 
of Fe ions of bcc Fe compared to hexagonal Gd2Fei7^. The same arguments support lower 
saturation magnetization value for Gd2Fei7. Thus even for highest Tc compound Gd2Fei7 
of i?2Fei7 series magnetic properties are found to be inferior to ones of bcc Fe. 



15 



V. ACKNOWLEDGMENTS 



We thank M.A. Korotin and V.I. Anisimov for helpful discussions. This work was sup- 
ported by the Russian Foundation for Basic Research, (project nos. 08-02-00021, 08-02- 
00712, and 09-02-00431), the Presidium of the Russian Academy of Sciences (RAS) "Quan- 
tum macrophysics" and of the Division of Physical Sciences of the RAS "Strongly correlated 
electrons in semiconductors, metals, superconductors and magnetic materials". A.L. ac- 
knowledges "Dynasty" Foundation. I.N. acknowledges Russian Science Support Foundation 
and MK-2242.2007.2. 



1 X. C. Kou, F. R. de Boer, R. Grossinger, G. Wiesinger, H. Suzuki, H. Kitazava, T. Takamasu, 
and G. Kido, J. Magn. Magn. Mater. 177-181, 1002 (1998). 

2 K. H. J. Buschow, in Handbook of Magnetic Materials (Amsterdam, Elsevier, 1991), Vol. 10, 
463. 

3 D. Givord and R. Lemaire, IEEE Trans. Magn. MAG-10, 109 (1974). 

4 H. Sun, J. M. D. Coey, Y. Otani, and D. P. F. Hurley, J. Phys.: Condens. Matter 2, 6465 
(1990). 

5 D. B. dc Mooij and K. H. J. Buschow, J. Less-Common. Met. 142, 349 (1988). 

6 T. H. Jacobs, K. H. J. Buschow, G. F. Zhou, X. Li, and F. R. de Boer, J. Magn. Magn. Mater. 
116, 220 (1992). 

7 B. G. Shen, F. W. Wang, L. S. Kong, and L. Cao, J. Phys.: Condens. Matter 5, L685 (1993). 

8 R. van Mens, J. Magn. Magn. Mater. 61, 24 (1986). 

9 Z. W. Li, X. Z. Zhow, and A. H. Morrish, Phys. Rev. B 51, 2891 (1995). 

A. G. Kuchin, A. S. Ermolenko, and V. I. Khrabrov, Phys. Met. Metallogr. 86, 276 (1998). 

1 B. Shen, B. Liang, Z. Cheng, H. Gong, S. Zhang, and L. Zhang, J. Phys. D: Appl. Phys. 31, 
2438 (1998). 

2 Yu. V. Knyazev, A. V. Lukoyanov, Yu. I. Kuz'min, A. G. Kuchin, and I. A. Nekrasov, Phys. 
Rev. B 73, 094410 (2006). 

3 V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, J. Phys.: Condens. Matter 9, 767 
(1997). 



16 



14 O. K. Andersen, Phys. Rev. B 12, 3060 (1975). 

15 Indeed, for the hexagonal phase of Gd2Fei7 there are two inequivalent types of Gd atoms. 
From calculating of Gd atomic sphere radii through merging of Hartree potentials one 
can obtain R(Gdl)=2.86 a.u. and R(Gd2)=3.89 a.u. which are quite different from each 
other. Since Gd atoms do not play leading role in our computation we chose for simplicity 
R(Gdl)=R(Gd2)=2.86 a.u. by introducing additional empty spheres. It does not affect signifi- 
cantly in its turn Fe atomic sphere radii which are more important here. 

16 K. H. J. Buschow, Rep. Prog. Phys. 40, 1179 (1977). 

17 O. Gunnarsson, O. K. Andersen, O. Jepsen, and J. Zaanen, Phys. Rev. B 39, 1708 (1989). 

18 V. I. Anisimov and O. Gunnarsson, Phys. Rev. B 43, 7570 (1991). 

19 A. I. Lichtenstein, M. I. Katsnelson, and G. Kotliar, Phys. Rev. Lett. 87, 067205 (2001). 

20 J. Hubbard, Phys. Rev. B 19, 2626 (1979). 

21 R. E. Prange and V. Korenman, Phys. Rev. B 19, 4691 (1979). 

22 C. S. Wang, R. E. Prange, and V. Korenman, Phys. Rev. B 25, 5766 (1982). 

23 J. Mathon, Phys. Rev. B 27, 1916 (1983). 

24 A. I. Liechtenstein, M. I. Katsnelson, V. P. Antropov, and V. A. Gubanov, J. Magn. Magn. 
Mater. 67, 65 (1984). 

25 V. V. Mazurenko, F. Mila, and V. I. Anisimov, Phys. Rev. B 73, 014418 (2006). 

26 J. P. Liu, F. R. de Boer, P. F. de Chatel, R. Coehoorn, and K. H. J. Buschow, J. Magn. Magn. 
Mater. 132, 159 (1994). 

27 Since distribution of atoms becomes more uniform with distance there is quite difficult to 
determine unambiguously upper border of second coordination layer. Thus for both phases we 
chose it to be 4.3 A. This value corresponds to the beginning of the second thin out spherical 
layer. 

28 R. F. Sabiryanov and S. S. Jaswal, Phys. Rev. B 57, 7767 (1998). 

29 S. V. Halilov, H. Eshrig, A. Y. Perlov, and P. M. Oppeneer, Phys. Rev. B 58, 293 (1998). 

30 M. Pajda, J. Kudrnovsky, I. Turek, V. Drchal, and P. Bruno, Phys. Rev. B 64, 174402 (2001). 

31 Averaged over first coordination sphere exchange interaction value for hexagonal Gd2Fei7 is 

computed as I = Jll(1)zil Sff4S(ijSi^{t) (4)244(1) = 99 - 6 K > where 4& and z ab are taken 
from Tab. [Ill 

32 Averaged over first coordination sphere Fe-Fe distance for hexagonal Gd2Fei7 is computed as 

17 



r;heX 



rii(l)zii(l)+ri 2 (4)z 12 (l)+...+r44(4)z44(l) 



2.574 A, where and z a b are taken from Tab. HH 



r 2ll(l)+Zl2(l) + --+244(l) 

33 Ch. Kittel, Introduction to Solid State Physics (Wiley, New York, 1966). 

34 To get saturation magnetization at T=0 of Fe ion sublattice in hexagonal Gd2Fei7 we took mag- 
netic moment per formula unit M(£Fe)=36.54 \ib and unitary cell volume $7=521. 4x 10~ 24 cm 3 



(Ref. 



12b. 



18 




Q Gd6c 

O Fe2 (9c/) 

O Fe4(18/7) 

© Fe3(18f) 
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FIG. 1: (Color online) Rhombohedral T/i2^ni7-type structure of GcbFeir. 
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FIG. 2: (Color online) Hexagonal T/^iVz^-type structure of Gc^Fe^. 
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